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ABSTRACT 

Chemical and photometric models of elliptical galaxies with infall of primordial gas, 
and subsequent ejection of processed material via galactic winds, are described. Ensur- 
ing that these models are consistent with the present-day colour-luminosity relation 
and the measured intracluster medium (ICM) abundances, we demonstrate that the 
initial mass function (IMF) must be significantly flatter (i.e. x ~ 0.80) than the canon- 
ical Salpeter slope (i.e. x w 1.35). Such a "top-heavy" IMF is in agreement with the 
earlier conclusions based upon closed-box models for elliptical galaxy evolution. On 
the other hand, the top-heavy IMF, in conjunction with these semi-analytic infall mod- 
els, predicts an ICM gas mass which exceeds that observed by up to a factor three, 
in contrast with the canonical closed-box models. Time and position-dependent IMF 
formalisms may prove to be a fruitful avenue for future research, but those presently 
available in the literature are shown to be inconsistent with several important obser- 
vational constraints. 
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1 INTRODUCTION 

Supernovae (SNe)-driven winds, and in particular their role 
, ' in setting the timescale for the cessation of bulk star forma- 
, tion, have long been recognised as important components 
5^ ' of elliptical galaxy formation/evolution models. Such wind 
models provide a natural framework in which the intrinsic 
elliptical galaxy colour-metallicity-luminosity (CML) rela- 
tionships can be established (e.g. Larson 1974; Arimoto & 
Yoshii 1987; Gibson 1996a, b). An inescapable consequence 
of such SNe-driven winds is the enrichment of the intergalac- 
tic medium with the products of SNe nucleosynthesis (e.g. 
Larson & Dinerstein 1975). 

Such wind "pollution" is particularly evident in the hot 
intracluster medium (ICM) of elliptical-rich galaxy clusters. 
Early, Type II SNe-dominated winds are the favoured mech- 
anism at play here, as evidenced by the strong cv-element en- 
hancement signatures in this hot gas (e.g. [0/Fe]> +0.2 - 
Mushotzky et al. 1996). The absolute mass of iron in a clus- 
ter's ICM can be shown to be related to the host cluster's 
integrated V-band luminosity tied up in ellipticals; i.e. 
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with loose groups and rich clusters populating the low and 
high luminosity extrema, respectively. While important,^ 
these ICM constraints do not, by themselves, provide much 
information on the underlying star formation or initial mass 
function (IMF) formalisms. 

Previous papers in this series (Matteucci & Gibson 
1995; Gibson & Matteucci 1997) have adopted the classic 
closed-box model for a galaxy's evolution and concluded that 
IMFs significantly flatter-than-Salpeter's (1955) canonical 
power-law slope x — 1.35 are required, in order to fit the 
ICM constraints noted above. Similar conclusions have been 
drawn by David, Forman & Jones (1991) Zepf & Silk (1996) 
and Loewenstein & Mushotzky (1996). 

In the past year though, interest has been piqued by 
the adoption of gas infall models for ellipticals (e.g. Tantalo 
et al. 1996; Kodama & Arimoto 1997), in lieu of the afore- 
mentioned closed-box scenarios. Such studies retain the SNe- 
driven wind framework but consciously restrict themselves 
to the spectrophotometric properties of their model galaxies, 
neglecting the consequences for enriching the ICM. Because 



where the mass and luminosity are in solar units (Arnaud 
1994). Similarly, a cluster's ICM gas mass 

M ICM 

can be 

related to its early-type galaxy V-band luminosity by 



For example, because of the a-element overabundance signa- 
ture, late-time Type la SNe-driven wind models are no longer a 
viable dominant ICM-pollution mechanism. 
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there exists a level of degeneracy in the present-day CML re- 
lationships, due to uncertainties in the IMF (e.g. Table 8 of 
Gibson 1997), star formation methodology, SNe energetics, 
etc, Tantalo et al. and Kodama & Arimoto used, partially 
for convenience, IMF power law slopes close to the Salpeter 
(1955) value - the former, x = 1.35; the latter, x = 1.20. 

The time is now ripe to assess the infall model in light 
of the heretofore neglected ICM constraints - in particular, 
the abundance ratios and absolute gas masses - in order to 
better appreciate its strengths and weaknesses. What fol- 
lows represents the first application of gas infall models to 
the ICM constraints, with answering the question - in the 
presence of gas infall, is the IMF in cluster ellipticals dom- 
inated by high-mass stars (i.e. top-heavy) , or does it follow 
the more canonical Salpeter distribution? - being the highest 
priority. 

To this end, in Section ^ we remind the reader of the 
primary motivation for considering the infall model in the 
first place. Following this, in Section ^, we present a grid of 
models which parallels that of Kodama & Arimoto (1997). 
We follow their numerical formalism, for convenience sake, 
although its similarity to that of Tantalo et al. 's (1996) 
means our conclusions are independent of this choice. As we 
will show, the favoured Kodama & Arimoto grid, with IMF 
slope x ~ 1.2, underproduces ICM iron by approximately 
a factor of five. We next present an alternate grid of infall 
models consistent with the ICM abundance constraints. As 
discussed in Section ^, consistent with the results of ear- 
lier closed-box models, only IMF slopes significantly flatter- 
then-Salpeter (i.e. x < 0.8) were capable of recovering the 
galaxy CML relations and the ICM abundance constraints, 
albeit at the expense of the ICM gas mass constraint. Our 
results are summarised in Section 0. 



2 WHY INFALL MODELS? 

Recall that our earlier papers (Matteucci & Gibson 1995; 
Gibson & Matteucci 1997) adopted the straightforward 
closed-box model for galaxy evolution. In particular, the star 
formation rate ip was assumed to vary directly with the gas 
mass, as 

V(t) = vM s (t), (3) 

with the timescale for star formation v varying in such 
a way that the present-day CML relations were properly 
established. In the closed-box model, star formation is a 
maximum at t = 0, steadily decreasing thereafter. Gibson 
& Matteucci (1997) showed that one can successfully re- 
cover the [(V-K),Lv] relation for ellipticals, and still honour 
all of the ICM abundance constraints, within the closed- 
box model framework, but only for IMF slopes significantly 
flatter-than-Salpeter's (1955) - e.g. a slope x « 1, as op- 
posed to Salpeter's canonical x = 1.35, was necessary. 

It has since become apparent (e.g. Tantalo et al. 1996; 
Kodama & Arimoto 1997) though that this same closed- 
box model for ellipticals does not appear to have the same 
success in recovering the ultraviolet (UV)-optical CML re- 
lations, a point not explored in the Gibson & Matteucci 
(1997) study. The reason for this failure is the overproduc- 
tion of low-metallicity (i.e. Z=0) stars during the initial in- 
tense star formation regime; the signature of such a com- 



ponent is not seen in the integrated spectra of old stellar 
populations (Worthey, Dorman & Jones 1996). 

Several alternate star formation formalisms, aimed 
specifically at avoiding this low-metallicity overproduction, 
have since been published - the previously noted Tantalo 
et al. (1996) and Kodama & Arimoto (1997) studies being 
the most noteworthy. These infall models are similar in spirit 
to those designed to avoid the analogous G-dwarf prob- 
lem encountered by solar neighbourhood-enrichment mod- 
els (e.g. Tinsley 1980), although for ellipticals, the infall 
timescales are, of course, substantially different to the slower 
accretion rates expected for spirals. Following Tantalo et al. 
and Kodama & Arimoto, we adopt an infall rate which par- 
allels the free-fall timescale. Star formation, in both models, 
is assumed to follow equation ^, but instead, with an ini- 
tial gas mass of zero, increasing in accordance with the as- 
sumed gas infall rate law. The accompanying rapid increase 
in global metallicity means that very few low-metallicity 
stars are formed, avoiding the G-dwarf problem. 

In summary, gas infall models for elliptical galaxies have 
been considered as viable alternatives to the classic closed- 
box model, primarily because they provide a simple solution 
to the overproduction of low metallicity stars which plagues 
the closed-box models - such a G-dwarf problem manifests 
itself most clearly in the closed-box model's problem in re- 
covering the UV-optical CML relations. 



3 ANALYSIS 
3.1 The models 

Using the photo-chemical evolution code of Gibson 
(1996a,b;1997), we constructed a grid of seven elliptical 
galaxy models consistent with the "metallicity sequence" of 
Kodama & Arimoto (1997, hereafter KA97). In particular, 
KA97 adopted an IMF slope mildly flatter (i.e. x = 1.20) 
than Salpeter's (1955) x — 1.35, with a mass range of 
0.1 — > 60.0 Mq. The cessation of star formation (i.e. at 
time iGw) was treated as a free parameter, and varied to 
ensure recovery of the present-day elliptical galaxy CML re- 
lations]^ Reasonable values for the parameters governing the 
timescales for star formation and gas infall were adopted (i.e. 
0.1 Gyr), independent of galaxy model. The adopted nucle- 
osynthetic yields are described in Gibson & Mould (1997). 

The first block of models in Table 1 represent our 
pseudo-KA97 grid (i.e. those labelled x — 1.20); unlike the 
data presented in their Table 2, to which the reader is re- 
ferred to for additional parameters not relevant to the dis- 
cussion at hand, we have included the masses of gas, oxy- 
gen, and iron ejected to the ICM, for each model (in units 
of 10 9 M Q , 10 6 M , and 10 6 Mq, respectively). The initial 

t It should be noted that KA97's models are systematically red- 
der in (V-K) than ours, by ~ 0.06 mag. Tracing the source of this 
discrepancy is not possible at this time, as their simple stellar 
population (SSP) colours are not currently available. Our grid is 
based upon the published Bertelli et al. (1994) isochrones, with 
the low mass extensions outlined in Gibson (1996a). For the al- 
ternate grid of models to be discussed shortly, we ensure self- 
consistency by maintaining this systematic 0.06 magnitude offset 
in (V-K). 
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Figure 1. Mass of iron ejected to the ICM as a function of 
present-day V-band luminosity. The models of Gibson & Mat- 
teucci (1997, GM97), Elbaz et al. (1995, EAV95), and Kodama 
& Arimoto (1997, KA97) are shown. The characteristic cluster 
galaxy V-band luminosity LT, is indicated. 



gas mass reservoir associated with each model (i.e. Mr, in 
units of 10 9 M©) is also provided in column 2, as it was 
not included in KA97's Table 2. Column 7 shows the resid- 
ual reservoir gas mass at the time of galactic gas expulsion 
tGW, again in units of 10 9 Mq. 

Of immediate concern is the mass of iron ejected from 
each of the model galaxies (column 6 of Table 1). In Figure 
|l|we show this mass of iron Mf c , for the Kodama & Arimoto 
(1997) models of Table 1, as a function of the present-day V- 
band luminosity (open squares). For comparison, the predic- 
tions of Gibson & Matteucci (1997, their "standard model") 
and Elbaz, Arnaud & Vangioni-Flam (1995) are also shown, 
with filled circles and crosses, respectively. We stress that the 
latter two models have been specifically designed to honour 
Arnaud's (1994) cluster V-band luminosity-ICM iron mass 

relation (equation ^)j^| It should be readily apparent from 
Figure ^ that Kodama & Arimoto's models systematically 
underproduce (by a factor of ~ 5) the necessary iron to re- 
cover Arnaud's relationship (equation [l]). 



1000 



100 




7.0 7.5 8.0 8.5 9.0 

log(t) [Gyr] 

Figure 2. Evolution of the star formation rate ip as a function 
of time t, for the seven "metallicity-sequence" models of Kodama 
& Arimoto (1997), with IMF slope x = 1.20 (solid curve). Also 
shown are the parallel grid of models with slope x = 0.80 (dot- 
ted curve), constructed a posteriori to recover Arnaud's (1994) 
cluster luminosity-ICM iron mass relation. The galactic mass of 
the respective models decreases from top-to-bottom, as per the 
corresponding entries in Table 1. 



models still underproduce iron£| The only ingredient which 
could be varied, and recover the Mp^ M -Lv relationship suc- 
cessfully, was the IMF sZopeJ^] 

Keeping all other input parameters the same as for the 
x = 1.20 grid of models, we varied the IMF slope x and 
the initial gas reservoir mass Mr, until a grid was gener- 
ated which honoured Arnaud's (1994) M|§ M - relation- 
ship (equation |l|) and retained consistency with the ellipti- 
cal galaxy [(V-K),Lv] relation. An IMF slope x = 0.80 was 
found to be necessary, in order to recover the mean of this re- 
lationship. A subset of these models are shown in the second 
block of Table 1. The galactic wind epochs for these models 
ranged from 120 to 65 Myrs (from most to least luminous). 
For comparison, Figure ^ shows how the star formation rate 
evolves for this subset (dotted curves), as well as those for 
the standard x — 1.20 models (solid curves) of Table 1. 



3.2 Infall and the IMF 

Our next step was to investigate mechanisms whereby Ar- 
naud's (1994) relationship could be honoured, but from 
within the infall framework. Varying the timescale for star 
formation, the gas infall rate, or the yield selection, does 
not alter the conclusion of the previous paragraph - i.e. the 

t While it is true that both Gibson & Matteucci's (1997) and 
Elbaz et al. 's (1995) models obey equation |l], only the former 
is consistent with the optical-infrared elliptical galaxy CML rela- 
tions (Gibson 1996a). 



3.3 Potential problems with the infall scenario 

Besides our conclusion regarding the necessarily steeper- 
than-Salpeter (1955) IMF slope of x w 0.80, there are two 

§ Invoking a later-time Type la SNe-driven wind component of 
such magnitude as to fully recover the ICM iron mass-cluster 
V-band luminosity relationship (equation |l|) , drives the predicted 
ICM [O/Fe] to ~ —0.5, a factor of ~ 5 — > lOx lower than observed 
(Mushotzky et al. 1996). 

^ Bimodal IMFs are also a possibility, but in their simplest form 
(Elbaz et al. 1995) do not recover the observed CML relations 
(Gibson 1996a). 
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[h P t] 

Table 1. Elliptical galaxy models constructed to conform to the metallicity sequence of Kodama & Arimoto (1997) - i.e. the block 
containing models with an IMF slope x = 1.20. The second block of models (i.e. IMF slope x = 0.80) was constructed a posteriori 
to ensure compatibility with Arnaud's (1994) cluster V-band luminosity-ICM iron mass relationship (equation 1). The masses of gas, 
oxygen, and iron ejected at time to the ICM (columns 4, 5, and 6) at time *gw are m units of 10 9 Mq, 10 6 Mq, and 10 6 Mq, respectively. 
Column 7 is the mass of gas remaining in the initial reservoir at time tGWi a l so m units of 10 9 Mq. See text for details. 
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other points which should be drawn from the model grid 

of Table l|] First, the standard x = 1.20 grid already lies 
consistently ~ 0.1 — > 0.3 dex below the mean observed 
[[< Z >]v,Mv] relation (e.g. see Figure 1 of Gibson & Mat- 
teucci 1997). The new x — 0.80 grid, which is now consis- 
tent with the ICM constraint of equation only exacerbates 
this shortcoming, being a further ~ 0.1 dex removed from 
the mean. Second, and perhaps more importantly, the flat- 
ter IMF grid requires an initial gas mass reservoir ~ 5x 
more massive than that for the steeper IMF grid, for galac- 
tic models of the same present-day luminosity. 

Following Section 7 of Elbaz et al. (1995), we then use 
the models of Table 1 to find that the predicted ICM gas 
mass to cluster luminosity is 



M' CM « 8L V , 



for the x = 1.20 grid, and 
A4 CM « 55L v , 



(4) 



(•») 



for the x — 0.80 grid. This predicted ICM gas mass includes 
the mass of gas ejected at tGW (column 4 of Table 1) and 
the residual, unincorporated, reservoir gas (column 7). An 
additional component from late-time (i.e. t > taw) winds 
or ram pressure stripping can be included but does not alter 
the coefficients shown in equations ^ and ^. The coefficients 
can be reduced by a factor of ~ 2 to recover the relation 
based upon the ejected gas mass at taw alone. 

For the x = 1.20 IMF, equation || implies that ~ 40 -> 
15% (loose groups to rich clusters) of the total ICM gas (i.e. 
equation ^) has been accounted for; ~ 60% of this has been 
processed and ejected from cluster ellipticals, and ~ 40% 
is associated with residual reservoir-gas. Q The remaining 



II The mass-to-blue light ratios (M/Lb) for the x = 0.80 models 
are ~ 10 — > 15% lower than the corresponding x = 1.20 grid. 
While interesting, the uncertainties in deriving the observational 
M/Lb-Lb relation make any further discussion regarding this 
difference moot. 

** These percentages have associated errors of ~ 20%, as we 
have fit simple power-laws to the gas mass— cluster luminosity 



~ 60 — > 85% (loose groups to rich clusters) of the observed 
ICM gas must therefore be assigned to some gas component 
which simply does not partake in the star formation process 
(either "actively" through galactic wind ejection, or "pas- 
sively", by being associated with a specific galactic halo's 
initial gas reservoir). This is perfectly in keeping with the 
results of Matteucci & Vettolani (1988), David et al. (1991), 
Elbaz et al. (1995), and Gibson & Matteucci (1997). 

On the other hand, equation [|, coupled with the x = 
0.80 grid, implies that ~ 300 — > 100% (again, loose groups 
to rich clusters) of the total ICM gas (i.e. equation El) has 
been accounted for! Again, this gas can be considered to 
be ~ 50% "processed-and-ejected" and ~ 50% "initial gas 
reservoir-associated" - no other cluster gas constituent need 
be invoked. The problem lies, of course, in that for any sys- 
tem "poorer" than the richest clusters, there is an accom- 
panying overproduction of gas - i.e. the grid of infall models 
constructed with the most conservative input parameters re- 
quire, and yield, more gas than is actually observed, for all 
but the richest clusters. Some of this discrepancy may be 
alleviated if the galactic winds are successful in overcoming 
not only an individual galaxy's potential well, but that of 
the group's as well, dispersing some fraction of the "excess" 
gas to the general intergalactic medium. 



3.4 Variable IMFs? 

Because of the closed-box model's failure to properly recover 
the UV-optical CML relations (recall Section B) and the in- 
fall model's problems with gas "overproduction" , it might be 
tempting, by process of elimination, to support the notion 
of time (and position)-variable IMFs, a la Padoan, Nord- 
lund & Jones (1997). Detailed models based upon Padoan 
et al. 's "universal" IMF formalism will be presented else- 
where, but we would be remiss if we did not at least temper 



relations in order to follow the procedure laid out in Section 7 
of Elbaz et al. (1995). Strictly, a numerical integration should be 
performed, as discussed already in Gibson & Matteucci (1997), 
but for simplicity's sake we adhere to the Elbaz et al. formalism. 
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such "default" support, by noting two of the weaknesses in 
the (currently) best available variable-IMF elliptical galaxy 
models. 

Chiosi & Bressan (1997) have presented their first (pre- 
liminary) grid of elliptical galaxy models, based upon the 
Padoan et al. (1997) variable-IMF formalism. These models 
were designed primarily to satisfy one projection of the fun- 
damental plane (specifically, M/Lb versus Mb), while still 
retaining consistency with the underlying galaxy CML re- 
lations and observed a-element stellar population overabun- 
dance. While successful in satisfying these constraints, there 
are two equally important observational constraints which 
are clearly violated. 

First, local dwarf ellipticals such as Fornax, NGC 185, 
and NGC 205, each host populations of planetary nebulae 
with [O/Fe] typically ~ +0.2, and ages presumably on the 
order of several Gyrs (Richer, McCall & Arimoto 1997). An 
inescapable result of the Chiosi & Bressan (1997) "canoni- 
cal" dwarf elliptical (dE) model (see their Figure 2) is that 
dE planetary nebulae, of this typical age, should have [O/Fe] 
of < —1.4, a factor of ~ 40 lower than that observed. 

Second, and perhaps more important, is the predicted 
ICM [O/Fe] from the Chiosi & Bressan (1997) grid. As their 
dwarf ellipticals never develop galactic winds, it is really 
only the ~L* galaxies which can contribute to the observed 
[0/Fe]> +0.2 (Mushotzky et al. 1996). Again, referring to 
their Figure 1, we can see that their ~L» ellipticals develop 
galactic winds after ~ 3.5 Gyrs, by which point the ISM 
[O/Fe] has been reduced to ~ — 1 (their Figure 2). This 
factor of ~ 10 — » 20 discrepancy between observed and the- 
oretical ICM [O/Fe] points to a significant deficiency in the 
existing elliptical galaxy variable-IMF models. 

Third, the Chiosi & Bressan (1997) model (again, 
coupling their Figure 1 with the adopted star formation for- 
malism) shows a star formation rate, in the core alone (where 
< 10% of the galaxy's luminous mass lies), ranging from 
i[> ~ 6000 -> 1500 M /yr from redshifts of ~ 5 to ~ 1. The 
total global star formation rate is ~ 2 — » 3x this core value. 
Such star formation rates at z ~ 1 (i.e. ip > 3000 Mq), re- 
gardless of IMF, are clearly at odds with the cosmological 
number counts (e.g. Chariot & Silk 1995), although a more 
formal analysis, in a cosmological context, is postponed for 
the time being. 

To end on a positive note though, the above arguments 
do not necessarily preclude the possibility that future, more 
sophisticated, variable-IMF models, may reconcile all of the 
galaxy CML relations and ICM constraints. It should be re- 
iterated that the Chiosi & Bressan (1997) results are still 
preliminary, and that the full parameter space still needs to 
be explored. What is clear though, is that the models pre- 
sented thus far are not necessarily the panacea they might 
at first appear to be. 



4 SUMMARY 

We have constructed a grid of gas infall models of elliptical 
galaxies, based upon the favoured sequence of Kodama & 
Arimoto (1997). For an IMF power-law slope x = 1.20, the 
models are entirely consistent with the observed, present- 
day, colour-luminosity relations for cluster ellipticals. While 
previous closed-box models were equally successful in recov- 



ering the optical-infrared relations, the infall models do have 
the advantage of a better fit to the ultraviolet-optical rela- 
tions (albeit, at the expense of a poorer fit to the metallicity- 
luminosity relation). We assess qualitatively the ability of 
Padoan et al. 's (1997) variable-IMF scenario, as adopted by 
Chiosi & Bressan (1997), to circumvent the closed-box/infall 
"problems" , and conclude that the extent models admittedly 
minimise some of these problems, but only at the expense of 
introducing new ones. 

For the first time, we have employed infall models as in- 
put to simple cluster enrichment scenarios. For the published 
grids, which appear to favour IMF slopes similar to that of 
the canonical Salpeter (1955) value (i.e. x « 1.2 — > 1.4), 
an inescapable conclusion is that the predicted mass of iron 
ejected to a cluster's ICM is a factor of ~ 5 below that 
observed. Only by substantially "flattening" the IMF (to 
x « 0.8), were we able to construct an "infall" grid which 
was consistent with both the galactic colour-luminosity re- 
lations and ICM iron mass-cluster luminosity relation, al- 
beit at (a) the expense of modestly worsening the galactic 
met allicity- luminosity relation, and (b) overproducing ICM 
gas for all but the richest clusters. This latter point may be 
indicative of a fundamental flaw in the infall picture, but we 
prefer to reserve judgment on this point until more sophis- 
ticated models have been explored. 

In the meantime though, we are able to conclude that, 
by adopting an already accepted (e.g. Tantalo et al. 1996; 
Kodama & Arimoto 1997) formalism for elliptical galaxy 
evolution based upon gas infall, the IMF must necessarily 
favour the formation of massive stars (i.e. be "top-heavy") if 
the models are to conform to the present-day CML relations, 
as well as the ICM abundance constraints. This agrees with 
our earlier studies based upon the more conventional closed- 
box model (Matteucci & Gibson 1995; Gibson & Matteucci 
1997). It is encouraging that two such disparate approaches 
(i.e. closed-box versus infall) to modeling the evolution of 
ellipticals both lead to the same conclusion regarding the 
veracity of the top-heavy IMF hypothesis. 
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